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Abstract 

Dispersive approach to quantum chromodynamics is applied to the assessment of hadronic contributions to elec¬ 
troweak observables. The employed approach merges the corresponding perturbative input with intrinsically nonper- 
turbative constraints, which originate in the respective kinematic restrictions. The evaluated hadronic contributions to 
the muon anomalous magnetic moment and to the shift of the electromagnetic fine structure constant at the scale of 
Z boson mass conform with recent assessments of these quantities. 
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Numerous strong interaction processes are governed 
by the hadronic vacuum polarization function Y\{q^). 
Its ultraviolet behavior can be studied within pertur¬ 
bation theory, whereas its infrared behavior is only 
accessible within various nonperturbative approaches, 
e.g., lattice simulations IMl, operator product expan¬ 
sion GHn, instanton liquid model lElEl, and others. 

Certain nonperturbative information about the low- 
energy hadron dynamics is contained within dispersion 
relations. The latter are widely employed in various is¬ 
sues of theoretical particle physics, for example, the pre¬ 
cise determination of parameters of resonances IIT4ll . the 
extension of applicability range of chiral perturbation 
theory Esma, the assessment of the hadronic light- 
by-light scattering EH, and many others (see, e.g.. 
Refs. EEHUD- 

The dispersion relations render the kinematic restric¬ 
tions on the relevant physical processes into the mathe¬ 
matical form and impose stringent intrinsically nonper¬ 
turbative constraints on the pertinent quantities. Among 
the latter are the function n(^^), which is defined as the 
scalar part of the hadronic vacuum polarization tensor 

= i f /xe'""(o|r{7^(x)4(0)}|o) = 

= i iqfiqy - gfivq^) n(?^)/(127rb, (1) 

related R(s) function 

R{s) - Im lim n(s + ie)/7r, (2) 
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which is identified with the R-ratio of electron-positron 
annihilation into hadrons, and the Adler function ll3^ 


£>(eb = - 


dn{-Q^) 
d\nQ^ ’ 


(3) 


with - —q^ > 0 and s - > 0 being the spacelike 

and timelike kinematic variables, respectively. 

The dispersive approach to QCD ll3Tl435l (its pre¬ 
liminary formulation was discussed in Refs. [371) 
merges the aforementioned nonperturbative constraints 
with corresponding perturbative input and provides the 
unified integral representations for the functions on 
hand: 


An(^^ q^) = qi) + 

{cr - q^ — q\\d(J 


i: 


p{(T) ln| 


a-- q^m^ - q^ / cr 


/^OO J 

R(s) - R®(i) + 8{s - m^) j p(cr)—^, 
Js cr 

cr — nP' da 


D{QP) = + 

Cy + mP 
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a + a 

Here SYl{q\ q^ = Uiq^) - n(q'g), nP = 4m^, 6{x) is 
the unit step-function [0(x) = 1 if x > 0 and 0(x) = 0 
otherwise], the leading-order terms read 


AnW(^2^ ql) = 2 


(p-iW(p ^ip()-iW(pQ 


tan^ ip 


tan^ ifQ 


= Q{s - irPfi - (OT^/i)j ^ , 
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= 1 -H 


3[l - ^1+r' sinh-'(^'/2^]r‘, (9) 













and pier) stands for the spectral density 


DPT 


1 d 

p(cr) = - —-Im lim p{o- - ie) 

71 diner £^o+ 


= Im lim d(-cr - ie)ln. 
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d r(cr) 
d In cr 
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In these equations sin^^ = Ini^, sin^(^o = 

^ Ini^i and p{q^), r{s), d(Q^) denote the strong cor¬ 
rections to the respective functions, see Refs. lEHia 
for the details. 

It is worth noting that the representations 
conform with the results of Bethe-Salpeter calcula¬ 
tions ||40]| as well as of lattice simulations gn . The 
Adler function ([^ agrees with its experimental predic¬ 
tion in the entire energy range ll^ l42l |4^ (the stud¬ 
ies of D(Q^) within other approaches can be found in, 
e.g.. Refs. IBHlQl). Additionally, the dispersive ap¬ 
proach has proved to be capable of describing OPAL 
(update 2012, Ref. (HI) and ALEPH (update 2014, 
Ref. 1521) experimental data on inclusive r lepton 
hadronic decay in vector and axial-vector channels in 
a self-consistent way IlMlISTl (see also Refs. (11155]). 

The perturbative part of the spectral density can be 
calculated as (see, e.g.. Refs. 11611121) 


I d 
dr^M) _ 

d\ncr 


Im lim p (a- - is) = 

e-)0+ ^ 

: - Im lim .,(-cr - ie), 

7T £-> 0 + 


( 11 ) 


that provides the respective perturbative input to the 
integral representations The latter are by con¬ 

struction consistent with aforementioned nonperturba- 
tive constraints and corresponding perturbative results 
and constitute the “dispersively improved perturbation 
theory” (DPT) expressions for the functions on hand. 
At the one-loop level Eq. (H) assumes a quite simple 
form, namely, Ppert(o-) = (4/y6o)[ln^(cr/A2) -H 
where y6o = 11 - 2nf/3, rif denotes the number of 
active flavors, and A is the QCD scale parameter. The 
explicit expressions for the spectral function (H) up 
to the four-loop level are given in Ref. Il5^ (recently 
calculated respective four-loop perturbative coefficient 
can be found in Ref. Il58l ). The perturbative spectral 
function (111 will be employed hereinafter. 

Note that in the massless limit (m = 0) for the case 
of perturbative spectral function Eqs. 0 and 0 
become identical to those of the “analytic perturbation 
theory” (APT) ifTSl (see also Refs. lfT9H3T1l '). However, 
as discussed in Refs. (33H35l|43[55l, the massless limit 
loses the substantial nonperturbative constraints, which 
relevant dispersion relations impose on the functions 



Figure 1: Hadronic vacuum polarization function within various ap¬ 
proaches: DPT [Eq. fTl) , solid curve], APT [Eq. jl3h dashed curve], 
PT [Eq. jl41, dot-dashed curve], and lattice data (Ref. circles). 


on hand, that appears to be essential for the studies of 
hadron dynamics at low energies. 

In what follows it is convenient to employ the sub¬ 
tracted at zero form of Eq. 0, specifically 

n(2^) = An(o, -Q^) = An™(o, -q^) + 



As one can infer from Fig. the obtained hadronic 
vacuum polarization function (solid curve) is in a good 
agreement with lattice data ll59l (circles) (the rescal¬ 
ing procedure described in Refs. 16011^ was applied). 
The presented result corresponds to the four-loop level, 
A = 419 MeV, and rif - 2. Figure [T] also displays 
the one-loop Eq. (0 in the massless limit, which cor¬ 
responds to APT (dashed curve) 


An«(-e^,-eo = in||^j + -in 








(13) 


and the one-loop perturbative approximation of Y[{q^) 
(dot-dashed curve) 


An«(-e^,-eo = in(-j.-in 






«pert(G') 


In these equations a(Q^) = a(Q^)/3Q/(4n), 
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(14) 


^ /Solnz’ " A2 
is the one-loop perturbative running coupling, and 
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Figure 2: The subtracted muon anomalous magnetic moment (Aa^ = 
- ao, ao = 11659 x 10“^): theoretical evaluations (circles) and 
experimental measurement (shaded band). 


stands for the one-loop infrared enhanced analytic run¬ 
ning coupling ||62l|63l, which was independently redis¬ 
covered in Refs. Il64ll65l . 

The perturbative approximation of n(q^) (14 1 con¬ 
tains infrared unphysical singularities, that makes it 
inapplicable at low energies. The expressions (12 1 
and ( [T3| ) contain no unphysical singularities, but their 
infrared behavior is quite different. Specifically, the 
APT prediction (13 i diverges at ^ 0 (that makes it 
also inapplicable at low energies), whereas the DPT ex¬ 
pression 0 vanishes in the infrared limit and proves 
to be applicable in the entire energy range. 


The persisting few standard deviations discrepancy 
between the experimental measurements ll^ iTOll and 
theoretical evaluations EUlTll of the muon anomalous 
magnetic moment - 2)/2 makes the latter a 

challenging issue of particle physics. The uncertainty 
of theoretical estimation of is largely dominated by 
the leading-order hadronic contribution ITT^ 



which involves the integration of Y[{q^) over the range 
inaccessible within perturbation theory. 

The DPT expression for Tl(q^) (|^ contains no un¬ 
physical singularities and enables one to perform the 
integration in Eq. ( [T7| i without invoking experimental 
data on R-ratio, that eventually results in llTSll 

= (696.1 +9.5) X 10^'“. (18) 
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Figure 3: Theoretical evaluations of hadronic contribution to the shift 
of electromagnetic fine structure constant at the scale of Z boson mass. 


Ref. ITMI . The obtained estimation ( [T8] l appears to be in 
a good agreement with its recent assessments Il66l468l . 

The complete muon anomalous magnetic moment 
comprises the QED contribution ifTSll . the electroweak 
contribution GU, as well as the higher-order ll^ 
and light-by-light Gil hadronic contributions, that, to¬ 
gether with lead to = (11659185.1 + 10.3) x 

10 see Ref. 1351. The discrepancy between this value 
and experimental measurement = (11659208.9 ± 
6.3) X 10 GUI 1781 corresponds to two standard devi¬ 
ations. As one can infer from Eig. S the obtained 
conforms with its recent evaluations 16614681 . 

Another observable of our interest is the electromag¬ 
netic running coupling 

= a[l - AQ',,p(^^) - \ (19) 

which plays a key role in a variety of issues of precision 
particle physics. The leptonic contribution to 

Eq. ( |T^ can be calculated by making use of perturbation 
theory l80l . whereas the hadronic contribution 


a 


AQ'had(? ) = 'P 


pf 

UnP- 


R(s) d s 


■ q^ s 


( 20 ) 


involves the integration over the low-energy range 
and constitutes the prevalent source of the uncertainty 
of ai.^{q^), see, e.g.. Refs. 166][81]. 

To evaluate the five-flavor hadronic contribution to 
the shift of the electromagnetic fine structure constant at 
the scale of Z boson mass within DPT we shall follow 
the very same lines as above, that eventually yields 1351 

^«had(^z7 = (274.9 + 2.2) X 10-4. (21) 


This equation corresponds to the four-loop level and the 
quoted error accounts for the uncertainties of the param- 


This equation corresponds to the four-loop level and the 
quoted error accounts for the uncertainties of the param¬ 


eters entering Eq. (17i, their values being taken from eters entering Eq. (20i, their values being taken from 
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Ref. 1741 . The obtained estimation of (21 

is in a good agreement with its recent evaluations [66 
IMIESI, see Fig.[^ At the same time, Eq. ( |2T] ) together 
with leptonic 1(801 and top quark Il82l contributions re¬ 
sults in - 128.962 + 0.030, that also conforms 

with recent assessments of this quantity Il^l68ll79l . see 
Ref. llTSl for the details. 

The author is grateful to G. Bali, R. Kaminski, 
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for the stimulating discussions and useful comments. 
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